Radiative heat transfer between two dielectric nanogratings in the scattering approach 
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We present a theoretical study of radiative heat transfer between dielectric nanogratings in the 
scattering approach. As a comparision with these exact results, we also evaluate the domain of 
validity of Derjaguin's Proximity Approximation (PA). We consider a system of two corrugated 
silica plates with various grating geometries, separation distances, and lateral displacement of the 
plates with respect to one another. Numerical computations show that while the PA is a good 
approximation for aligned gratings, it cannot be used when the gratings are laterally displaced. We 
illustrate this by a thermal modulator device for nanosystems based on such a displacement. 
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Recent experiments and theoretical work have given 
promising perspectives in the field of radiative heat trans- 
fer in the micrometer range [1-3 . It has been shown 
that radiative heat transfer greatly exceeds the black 
body limit for distances shorter than the average thermal 
wavelength, which is understood as an effect arising from 
the contribution of the evanescent waves. The studies of 
near-field heat transfer are of great interest to the design 
of both NEM S and MEMS which are naturally affected 
by possible side-effects of heat exchange at the nanoscale. 
Other potential applications lie in the fields of nanotech- 
nology, photonic crystals [4], metamaterials [SJE], ther- 
malphotovoltaics [7] [8] , multilayered structures [9] , im- 
proved resolution in nano-structure imaging, and new 
nano-fabrication techniques. 

While radiative heat transfer beyond Stefan- 
Boltzmann's law was observed experimentally [10] 
and described theoretically [11] over forty years ago, 
radiative heat transfer between two parallel flat plates at 
the nanoscale has been considered experimentally only 
recently [T2HT4] . The most interesting features of this 
field are the possible side-effects of non-trivial geometries 
on the thermal emission of nano-objects. Thus an in- 
depth study of heat transfer for different configurations 
has been performed over the years, ranging from the 
case of a particle facing a surface [T5HT8] , to particles or 
nanospheres facing each other [T9H23] . or more recently 
to the sphere-plane geometry [2] [24]. One should also 
note that for nearly fiat surfaces where roughness is 
considered as a perturbation factor, certain perturbative 
approaches can be used [25H27] . But for larger geo- 
metrical irregularities, more accurate methods become 
necessary [28] . These more complex geometries are best 
described through a scattering approach [24] I29H5T] . 
Another exciting perspective is the study of the variation 
in heat transfer brought forth by surface polaritons in 
certain materials [18] . In this paper we focus on the 
interplay between the surface waves excitation and the 
surface profile, as shown in Fig.l. 

The fact that the radiative heat transfer in near- 



field considerably changes with variation of the separa- 
tion distance between plane surfaces has already been 
shown [TJ [2] [32] . When introducing a profile for the in- 
terfaces, the flux is expected to depend on the relative 
lateral displacement of the two surfaces denoted 5, as 
seen in Fig.l. This is all the most interesting as a simple 
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FIG. 1: Two identical gratings facing each other at a distance 
L and relatively shifted by a lateral displacement 5. The 
corrugations have a period d, height a and thickness p . The 
filling factor p = p /d is given as a percentage of the period 
d. 

argument based on the Proximity Approximation sug- 
gests a strong modulation of the flux. Indeed, by assum- 
ing that one can use locally the plane-plane heat transfer 
coefficient, it is seen that the flux is maximum for S = 0. 
The validity of the proximity approximation has been 
discussed in the context of a plane-sphere [2J and be- 
tween two spheres [19] [20] . While it appears to be valid 
for radius of curvatures larger than 20/im, this validity 
in the context of lamellar gratings with sub- wavelength 
periods remains an open question. Here, we investigate 
this issue by using the exact formalism of scattering the- 
ory. Furthermore, we discuss the physical phenomena 
involved and show that the nature of the material needs 
to be taken into account when discussing the validity of 
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the Proximity Approximation. 

Based on the scattering formalism developed in [29] , we 
consider two corrugated profiles at temperatures T\ and 
T 2 , as shown in Fig.l. The heat transfer is constructed 
from the statistical average of the (x, y) sum over the z 
component of the Poynting vector S z and is thus related 
to a flux. We define the wave vector k = (k_L, k z ) with 
k z = \J~uP 1 /c 2 — k^_ defined with — tt/2 < arg/c 2 ^ n/2. 

Following [33], we then introduce the reflection oper- 
ators Ri(cj) and R 2 (cj) of the two gratings separated 
by a distance L, by which we understand the distance 
of closest approach, equal to zero at contact. We then 
set our scattering operators such that Si = Ri(w) and 
S 2 = e lkzL Il2(oo)e zkzL . According to the scattering for- 
malism for gratings developed in [33, 34], the scattering 
matrices are of dimensions 2 (27V + 1), where TV is the 
order of diffraction. 

We now define the operators J^ ,eui = \k n z W"l e " as 
constructed from the projectors on the propagative and 
evanescent sectors, respectively : 



U puj - 5 

u eu , — s 
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- sgn(cj 2 /c 2 
■ sgn(cj 2 /c 2 



(i) 

(2) 



where a = s,p represents the transverse electric and 
transverse magnetic polarizations, respectively. The 
thermal energy density per field mode at temperature 
T writes ct(oS) = Huu/ { e huJ / kBT — l). We can now ex- 
press the heat transfer coefficient between two gratings 
of same corrugation depths a as : 



1 



/ 



duo 
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(e Tl (cu) - e T2 (uj))H 12 
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with 



+7r/d 



4tt 2 



tr(DWiD t W 2 ) (4) 
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k x — — 7r/d J /c y Gff 

D = (1-S 1 S 2 )- 1 (5) 
= Er i -S 1 Er i S 1 t + S 1 5]!. a; 1 -5]!. a iS 1 t (6) 
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Ss^f S 2 + Sa^f - £f S 2 (7) 



It is noteworthy that the heat transfer depends on 
the shape and material properties of the gratings only 
through their scattering matrices Si and S 2 . Further- 
more, the factor 6^(00) — eT 2 (oj) introduces a cut-off for 
all frequencies larger than fc^T/fi. It is hence H\ 2 in 
equation Q , which corresponds to the sum of the trans- 
mission factors of the modes, that gives rise to the inter- 
esting modes pertaining to the near-field contribution. 

Note also that the first perpendicular wave vector com- 
ponent k x belongs to the first Brillouin zone between 
—ir/d and +7r/d, whereas k y E M is not restricted. A 
practical challenge of the numerical integration of h lays 
in the choice of the boundaries of uo and k y through a 
careful study of the integrand of equation Q plotted over 
the whole range of frequencies to determine the modes. 
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FIG. 2: Heat transfer coefficients as a function of separation 
distance L, when the gratings are not laterally displaced (blue 
solid line) and when they are by half a period (blue dashed 
line). This is compared with the Proximity Approximation 
in red. The gratings have a period d = 1500nm, filling factor 
p = 20%, and groove depth a = 500nm. 
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FIG. 3: Field modulus map of a given source dipole placed 
in the middle of a corrugation and right under the surface. 
The field is here represented only in the upper grating, so 
as to highlight where the absorption takes place. This is for 
gratings at a separation distance L = 25nm. The two figures 
on the left display the profiles in the xz— plane (in green) when 
they are aligned (6 = 0), and the two figures on the right when 
they are laterally displaced by half-a-period (5 = d/2), both 
for two different wavelengths A = 8.75/xm (top) and 9.15/im 
(down). 



We will from now on consider two gratings of silica 
glass Si02, the dielectric properties of which are given 
in [35]. This material is chosen as it supports surface 
phono n-polaritons, which are known to enhance the flux. 
The gratings temperatures are supposed to be T\ = 310 
K and T 2 = 290 K. Two sets of data are systematically 
computed : the first corresponds to zero lateral displace- 
ment of the two plates along the x-axis (5 = 0) so that 
the corrugation maxima directly face those from the op- 
posite profile. The second corresponds to a lateral dis- 
placement of half the grating period (5 = d/2), so that 
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the corrugations maxima face the corrugation trenches 
of the opposite profile. In near-field, the two structured 
plates expose a larger surface to each other at 5 = 
than at 6 = d/2, so that we expect a strong modulation 
of the heat transfer coefficient which will be discussed 
later. This is based on the assumption that the plane- 
plane heat transfer coefficient is locally valid. 

The results of the scattering approach can be compared 
with the PA, which consists on the weighted sum of the 
planar normal contributions ho(L) depending on the local 
separation distances L within each period. Assuming 
that p < 50%, we have for S <p' : 

h* s K {L)= ^^h (L) + ^-ho(L + a) (8) 
+ h (L + 2a) 

For S > p, we find the following saturation value of : 

h P /(L)= *£ho(L + a) (9) 

In what follows, we study in detail the interplay be- 
tween surface waves and corrugations. In particular, we 
assess the validity of the PA. Fig. 2 shows the heat trans- 
fer coefficient for 6 = and 6 = d/2, as a function of 
the separation distance L, for two gratings of period 
d = 1500nm, filling factor p = 20%, and groove depth 
a = 500nm. The results are compared with the PA. Re- 
gardless of the distance, we can see that the PA is a good 
approximation to the heat transfer coefficient at S = 0, 
but not at S = d/2. At L = 25nm, the error of the PA is 
of - 3% for 5 = 0, and of - 35% for S = d/2. 

The reason for this is illustrated on Fig. 3, which shows 
the field modulus map for a given source dipole that is 
placed in the middle of a corrugation right under the 
surface, and which is oriented perpendicular to it. The 
color scale is logarithmic. The intensity or square mod- 
ulus of the electric field is represented only in the upper 
grating so as to highlight the place of absorption. The 
gratings have a separation distance L = 25nm, corruga- 
tion depth a = 500nm, period d = 1500nm, and filling 
factor p = 20%. Two different wavelengths A = 8.75/im 
and A = 9.15/im are considered, knowing that Si02 has 
two resonance frequencies at A = 8.75/im and A = 21 /am. 
In the case where 5 = and A = 8.75/im, we see that 
the field is clearly both intense and confined. As 8.75/im 
corresponds to the horizontal asymptote of the surface 
phonon dispersion relation, a large number of modes with 
different values of the wave vector are excited. This leads 
to a highly localized subwavelength hot spot. At 9.15/im, 
the spot is broader as expected : this is similar to the loss 
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FIG. 4: Heat transfer coefficients as a function of the separa- 
tion distance L between two plane mirrors of SiCb (red solid 
curve), compared with the black body limit (blue dashed line). 
One can divide the separation distance in three domains A, 
B, and C, respectively corresponding to the extreme near- field 
below 200nm, to the near-field from 200nm to 10/im, and to 
the domain of Stefan-Boltzmann's law beyond 10/xm. This 
can be seen by the change of the slope of the curve along 
these three ranges. 

of resolution of surperlens away from the resonance. On 
the right column of the figure, we show the intensity for 
S = d/2. It is seen that the heated region is delocalized 
so that PA is clearly not valid. In this regime, the heat 
transfer is no longer due to a dipole-dipole interaction 
through the gap. Instead, a dipole excites modes of the 
structures. In turn, these spatially extended modes pro- 
duce dissipation in the walls. This discussion indicates 
that PA is valid if the gap width does not vary signifi- 
cantly on a length scale given by the spatial extension of 
the modes. 

To further illustrate this qualitative dependence of the 
radiative heat transfer on separation distance, we show 
in Fig. 4 the heat transfer coefficients as a function of 
the separation distance L between two plates of Si02- 
One can distinguish three domains A, B, and C, cor- 
responding respectively to the extreme near-field below 
200nm, to the near-field from 200nm to 10/im, and to the 
domain of Stefan-Boltzmann's law beyond 10/im. The 
heat transfer coefficient changes in slope along these three 
ranges: the strongest contributions come respectively 
from the dipole-dipole interaction, from surface phonon- 
polaritons, and from the classical radiative heat transfer. 
The contribution in the first domain corresponds to the 
localized heat transfer seen in the upper-left map of Fig. 3, 
whereas the main contribution in the second domain cor- 
responds to the delocalized heat transfer mediated by the 
surface wave seen on the right maps of Fig. 3. 

It is also instructive to study the heat transfer modu- 
lation as a function of the corrugation period d, as shown 
in Fig. 5. We have selected six types of gratings with cor- 
rugation periods ranging from d = 250 to 1500nm, each 
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FIG. 5: Heat transfer coefficients as a function of grating pe- 
riod d, when the gratings are not laterally displaced (solid blue 
line) and when they are displaced by half a period (dashed 
blue line). This is compared with the PA in red. The gratings 
have a groove depth a = 500nm, filling factor p = 20%, and 
are at a separation distance L = lOOnm. 
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FIG. 6: Heat transfer coefficients as a function of fill- 
ing factor p, when the gratings are not laterally displaced 
(solid blue line) and when they are displaced by half a pe- 
riod (dashed blue line). Respective PA predictions are in red. 
The dotted gray line is the percentage of the modulation fac- 
tor hs=o/hs=d/2- Gratings have a period and groove depth of 
500nm, and are separated by a distance L = 100 nm. 



with a groove depth a = 500nm and filling factor still 
fixed at p = 20%. The separation distance is L = lOOnm. 
The fact that the heat transfer coefficients at 5 = does 
not vary much with a change of period is a further con- 
firmation of the validity of the PA in this configuration. 
At 5 = d/2, however the scattering and PA results radi- 
cally differ for small periods, but tend to agree for large 
periods. The reason for this is that when d — >• oo, the 
ratio a/d tends to zero and we expect the heat transfer to 
be well approximated by the plane-plane case, and hence 
the PA. 

Let us finally turn to the discussion of the modulation 
effect. Fig. 2 shows that the heat transfer depends dra- 
matically on the lateral displacement of the two surfaces, 
opening the possibility of a strong modulation via only 



lateral displacement of one of the two plates at a fixed 
distance. To assess the possible performance of such a 
system as a thermal modulator ??, we investigate the 
modulation factor /^=o/^5=d/2 f° r different filling fac- 
tors. The results are illustrated in Fig. 6 for gratings with 
a period and groove depth a = 500nm, and a separation 
distance L = 100 nm. For these large separations, the 
modulation factor hs=o/h§ =d / 2 still reaches a maximum 
of about 2.2, at a filling factor corresponding to 20% of 
the total grating period. At short distances (L ~ 25nm) 
it can reach up to 35 (c.f. Fig. 2.). 

So based on a - to our knowledge first - study of ra- 
diative heat transfer for corrugated dielectric plates - we 
have clarified the origin of the success and failure of the 
PA by analyzing the interplay between surface waves 
resonances and corrugations, we have shown for vari- 
ous nanograting geometries and separation distances that 
Derjaguin's Proximity Approximation is clearly valid for 
6 = but fails for 6 = d/2. The key to the understand- 
ing of the system is the comparison of the lateral length 
scale of the surface corrugation with the lateral extension 
of surface waves involved in the heat transfer. Finally, we 
have narrowed down the optimum geometrical parame- 
ters of a thermal modulator device for nanosystems based 
on a lateral displacement of two corrugated plates fac- 
ing each other at fixed distance. We found in general a 
stronger modulation for small filling factors and separa- 
tion distances, and for large grating periods. In certain 
regimes it is possible to reach a modulation factor of more 
than 35. An in-depth study of the modes accounting for 
the most important part of the heat transfer would be 
an interesting prospect as well as to further enhance the 
modulation by using a broader range of materials [36] 
such as different alloys combining the polaritons of cer- 
tain dielectrics and the near-field properties of metals. 
The issue of heat transfer in near-field in the case of 
coatings [9], phase change materials [32j [36], metama- 
terials [6] [37] , or graphene-covered dielectrics [3] in this 
regard should also be explored. 
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